NOT TO BE TAKEN FROM THIS ROOM 


id 


Sat Pal ECR OES 


FeO ad hones Sune ns te ten 


For Reference 


NOT TO BE TAKEN FROM THIS ROOM 


Digitized by the Internet Archive 
in 2018 with funding trom 
University of Alberta Libraries 


https://archive.org/details/rheologyofsuspben00aema 


Jigen 
19 a F) 


452, 
THE UNIVERSITY OF ALBERTA 


THE RHEOLOGY OF SUSPENSIONS OF SPHERES 


EFFECT OF SIZE DISTRIBUTION 


A THESIS 
SUBMITTED TO THE FACULTY OF GRADUATE STUDIES 
IN PARTIAL FULFIIMENT OF THE REQUIREMENTS FOR THE DEGREE OF 


MASTER OF SCIENCE IN CHEMICAL ENGINEERING 


DEPARTMENT OF CHEMICAL AND PETROLEUM ENGINEERING 
by 


A.E. MATHER, B.Sc. 


EDMONTON , ALBERTA 


OCTOBER, 1962 


ABSTRACT 


The effect of the size distribution of the disperse phase 
upon the rheological behavior of suspensions of spheres was studied 
with a rotational viscometer. The suspensions were found to behave 
as Newtonian fluids over the range of concentration studied (0 to 50 
volume percent). For the three unimodal suspensions studied (average 
particle diameter 47, 86.5 and 390 microns), no effect of particle 
size upon the viscosity was found. A bimodal distribution of part- 
icle sizes was found to have a significant effect upon the viscosity 
of the suspensions, the viscosity being lower than that of a unimodal 
suspension at the same concentration. It was found that the minimum 
viscosity for any volume concentration occurs when 70% of the total 
solids are large particles and 30% are small particles. Increasing 
the diameter ratio of large to small particles also decreases the 
viscosity at any volume concentration. These facts are explained 
in terms of the packing of the smaller spheres into the interstices 
between the larger spheres. The results of this work agree well 


with those of Sweeny (46). 
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I. INTRODUCTION 


This work is an experimental investigation of the effect of 
the size distribution of the solid phase upon the rheological behavior 
of suspensions of spherical particles. The ultimate aim of work such 
as this is the thorough understanding of the process of flow of part- 


iculate solids in liquids. 


The flow behavior of solid-liquid mixtures is important in 
a number of industrial processes, especially those involved in the 
transport of raw materials such as cement, sand, gravel, coal and 
various ores. The sucessful application of suspension transport 
depends to a large extent on the ability to predict the power consum- 
ption involved in the process. This in turn is directly dependent 
upon the rheological behavior of the solid~liquid mixture. It would 
be desirable to be able to relate the physical properties of the 
system - the particle size, shape, size distribution and concentration 
of the disperse phase, to the rheological behavior of the suspension. 
At present, it is not possible to do this because of the scarcity of 
both theoretical and experimental work on the effect of the various 


physical properties upon the rheological behavior of the suspension. 


A survey of the literature showed that even suspensions 
of spheres, the most idealized system possible, could not be fully 


characterized by the physical parameters. The effect of the size 
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distribution of the disperse phase had not been studied extensively, 
although there was some indication that a wide size distribution 
contributed to a lower consistency than a narrow size distribution 

at the same concentration of solid phase. At the time that this work 
was begun, there were only a few cases of study of the effect of size 
distribution upon the rheological properties of suspensions of spheres 
and these were at low concentrations. For this reason, it was 

decided to undertake an experimental investigation of the effect of 
the size distribution of the spheres upon the rieological behavior of 


the resulting suspension. 


The rheological properties of the suspensions used in this 
investigation were measured by means of a concentric cylinder visco~ 
meter. This type of viscometer was chosen because of its advantages 
over other instruments in the measurement of rheological properties 
of suspensions. These include the nearly constant shear rate through- 
out the sample, the ease of operation and the ability to detect and 


study thixotropic breakdown in a fluid. 


This study to determine the effect of the size distribution 
upon the rheological behavior of suspensions of spheres is presented 
under the subjects of rheological and viscometric theory, a literature 
review of the most important work on suspensions of spheres, in 
particular that concerning the effect of size distribution, the experi= 
mental equipment and method used in this investigation, experimental 
results and a discussion and some conclusions warranted as a result of 


this work. 
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II. THEORY AND LITERATURE REVIEW 


1. Rheological Background 


Fluid Classification 


Fluids are classified as Newtonian or non-Newtonian depending 
upon their behavior, at constant temperature and pressure, under imposed 
shearing stresses which give rise to laminar motion. A Newtonian fluid 


is one which obeys Newton's (36) law, 
Sigmar 3 cre & 


where a is the coefficient of viscosity, often simply termed 


viscosity, poise (gm/sec-cm) 


ae 2 
| is the unit shearing stress, dynes/cm 

du : Sil 

ax is the shear rate, sec 


A graphical representation of the flow curve of a Newtonian fluid is 
shown in Figures 1 and 2 on arithmetic and logarithmic coordinates. 
The flow curve of a Newtonian fluid is completely characterized by 


the single physical property of viscosity. 


By exclusion, fluids which do not obey Newton's law are 
called non-Newtonian fluids. Non-Newtonian fluids may be further 


subdivided into three broad categories: 
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i) Fluids with properties independent of time or duration 
of shear 
ii) Fluids with properties dependent upon the duration 
of shear as well as the shear stress. 

iii) Fluids with properties dependent upon the extent of 
deformation (strain) as well as shear stress. This 
last case involves fluids which have characteristics 
of solids and show elastic recovery. This category 


will not be considered further. 


Non-Newtonian fluids in categories i) and ii) are subdiv- 


ided into five general types: 


i) Bingham plastic ii) Thixotropic 
Pseudoplastic Rheopectic 
Dilatant 


These will now be considered in turn. 


Bingham plastic 


Fluids exhibiting Binghém plastic behavior do not flow until 
the shearing stress exceeds a minimum value. This minimum shearing 
stress, Ty » the yield stress, is the intercept on the shear stress 


axis in Figure 1. The equation of the flow curve is? 
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where Mr is the coefficient of rigidity and is the slope of the flow 
curve. On logarithmic coordinates, Figure 2, the flow curve becomes 
asymtotic to Ty at low shear rates and approaches a slope of unity 


at high shear rates. 


Bingham plastic behavior has been ascribed to the presence 
of a three-dimensional structure within the fluid while at rest, 
which collapses completely when the yield stress is exceeded. Actually, 
the Bingham plastic model represents a very idealized behavior and it 


is doubtful that any true Bingham fluids exist. 


Pseudoplastic 


From the flow curves of Figure 1, it is seen that the ratio 
of the unit shearing stress and the rate of shear for a non-Newtonian 
fluid depends upon the rate of shear. This ratio is different at 
different rates of shear and hence the consistency varies. To 
refer to these different consistencies as viscosities is incorrect; 
however, they may be referred to as apparent viscosities which have 


meaning only when the rate of shear is given. 


Many fluids, known as pseudoplastics, exhibit a type of flow 
in which the rate of flow increases faster than the shearing stress. 


Thus the apparent viscosity decreases as the rate of shear increases. 
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Pseudoplastic fluids exhibit no yield value and the flow curve becomes 


almost linear at very high rates of shear. 


Many explanations have been proposed to account for this 
behavior. One explanation is that highly asymmetric molecules or 
particles are randomly dispersed and entangled when the fluid is 
at rest. Under shear they tend to align themselves in the direction 
of the shearing motion. A second possible explanation of pseudoplas~ 
ticity may be the progressive shearing off of liquid carried by part- 
icles as envelopes of solvation, thus decreasing their effective size 


with increasing rate of shear. 


The logarithmic plot of shear stress versus rate of shear 
for these materials frequently approaches a straight line with 4 
slope between zero and unity. For this reason, an empirical rela= 
tion known as the power law is used to characterize fluids of this 


type. It is of the form: 
du \" r 
a x (gz) (3) 


This relation was originally proposed by Ostwald (37) and 


has been critically discussed by Reiner (39). 
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Dilatant 


Dilatant fluids are those for which the apparent viscosity 
increaseg with increasing rate of shear. Originally, the term dilat~ 
ant was applied to those systems which increased in volume, or dilated 
under shear. This type of behavior was found in concentrated suspen~ 
sions by Reynolds (40). When the suspension is at rest the voidage is 
a minimum 4nd the liquid is sufficient only to fill the voids. On 
shearing, the dense packing of the system is broken and the material 
"dilates'", increasing the voidage. There is now insufficient liquid 
to lubricate the flow of particles past each other and the applied 
stresses have to be much greater. This causes the apparent viscosity 
to increase rapidly with increasing rates of shear and is termed 


volumetric dilatancy. 


Freundlich and Roder, (17) in their work on starch~water 
and quartz-water suspensions, found that the apparent viscosity 
increased with increased shear rate. This was called dilatant 
behavior and the term dilatancy has been used to describe both 
volumetric dilatancy and rheological dilatancy. The latter can 
eecur at concentrations much below that required for volumetric 


dilatancy. 


Some dilatant systems cam be characterized rheologically 
by means of the "power law", Equation (3), except that in this case 


the exponent is greater than unity. 
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Thixotropic 


The classification, thixotropic, is used to denote those 
fluids whose consistency decreases with time under shear and builds up 
again when at rest. Thixotropy was first defined as an isothermal, 
reversible sol=gel~sol transformation, This definition has been 
broadened to include fluids that exhibit an isothermal reversible 
transformation from a higher to a lower consistency gel. A physical 
model for thixotropic fluids consists of coiled chains of long molecules, 
entangled in the fluid at rest. Upon shear, the chains begin to align 
themselves in the direction of the shear. At the same time breaking 
and disentanglement of some of the chains may occur, These processes 
are not instantaneous and the shear stress, at constant shear rate, is 


observed to decrease with time. 


Rheopectic 


Rheopectic fluids exhibit the opposite behavior to thixotropic, 
that is, at constant rate of shear, the consistency increases with time. 
Freundlich and Juliusberger (18) first used the term to describe the 
"solidification" of vanadium oxide suspensions under low rates of shear. 
It has been postulated that rheopexy is simply the same as dilatency 
except that a finite time interval is required. Rheopexy has been noted 


in various sols and in some bentonitic clay suspemsions (57). 
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2+ Theory of the Concentric Cylinder Viscometer 


The concentric eylinder or Couette viscometer is well suited 
for the study of the rheological properties of fluids, Its simple 
geometry allows the direct relation of the observed variables with 
the shearing stress and rate of shear. The shear rate is nearly 
constant throughout the sample and the instrument is simple to operate. 
This type of viscometer is also one of the few that can be used for a 
quéantitetive study of thixotropic breakdown because the shearing stress 
can be measured as a function of time at constant shear rate and 
temperature, Another reason for the suitability of Couette viscometers 
for measurement of the rheological properties of suspensions is that 
no wall effect is noted. Workers studying the rheology of suspensions 
of spheres in capillary viscometers have found that the viscosity of 
Newtonian suspensions is markedly dependent upon the capillary diameter 
(51), (24), (53). However, Eveson, Whitmore and Ward (14) found no 
effect of the gap width in a Couette viscometer upon the measured 
viscosity. For this reason, coaxial cylinder viscometers are to be 
preferred to capillary tube viscometers for the measurement of the 


rheological behavior of suspersicnss 


The Couette viscometer consists of two concentric cylinders 
arranged so that the liquid sample is placed between them. One 
cylinder, usually the outer, is rotated at a constant speed and the 


torque transmitted by the liguid to the immer cylinder is measured 
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by means of the deflection of a torsion wire.e The apparatus is usually 


placed in a thermostat to insure a constant temperature of the sample. 


It is possible to use the alternate system of the inner 
cylinder rotating and the outer cylinder fixed. Taylor (48) has 
shown that the outer cylinder rotating is considerably more stable 
than the inner cylinder rotating, in that laminar flow occurs up 


to much higher speeds. 


The experimental measurements consist of the angular velocity, 
f2., of the outer cylinder and the torque, M, exerted on the inner cyl- 
inder. From the measurement of these relations it is possible to 
determine the flow curve of the liquid. The following assumptions 
are made regarding the motion of the liquid: 

i) the liquid is incompressible, 

ii) the motion of the liquid is laminar, 

iii) the laminae are circles on the horizontal planes 
perpendicular to the axis of rotation, 

iv) the motion is constant with time, (there is no acceler- 
ation) 

v) there is no slip at the walls, ise. the velocity of the 

fluid at the surface of the cylinders is zero, and, 


vi) the cylinders are assumed to be infinite in length. 
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With these assumptions, it is possible to obtain the following 


equation (see Appendix G) which relates f& and shearing stress for any 


arbitrary time-independent fluid. 


where the rheological behavior of the fluid is expressed as 


Oe SVs 


and shearing stress is related to torque by 


Tos ae 


(4) 


(5) 


(6) 


where r = radial distance (some value between the radius of the inner 


cylinder and that of the outer cylinder) 


Dp 
Wt 


length of cylinder immersed in the sample. 


Equation (4) can be integrated if a £(T) is assumed for a 


fluid. For a Newtonian fluid, £(T) = T/m and 
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Expressing ay in terms of M: 
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where R; = radius of the inner cylinder 


I 


Ro = radius of the outer cylinder 


This equation is the well-known Margules equation. 


When dealing with non-Newtonian fluids, it is desirable to 


be able to construct the flow curve (shearing stress - rate of shear 


diagram) without the assumption of any particular behavior for the fluid. 


This is easily done for the shearing stress since Equation (6) is valid 
for any fluid. At the surface of the bob, therefore, the shearing 


stress for any fluid is: 


Te; = M (9) 
2TTRy bh 

To obtain the rate of shear at the surface of the bob for any fluid is 

more difficult. Pawlowski (38) and Krieger and Elrod (28) have solved 

this problem. The equations are independent of the rheological nature 

of the fluid. The Krieger and Elrod equation consists of an infinite 


series which may be truncated to three terms if the ratio s = Ro/Ri 


is less than 1.2. 
os = _4TN tot tea eo Wl ( L/n" ~ 1)7 + an" = We eel \ (10) 
=a i Uf f d log M 


where kj and kz are instrument constants 


ky = s*=-1 (1+ 2/3 (ins) ) (11) 
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n'' is the slope of a logarithmic plot of torque M versus rotational 
speed N at the value of N in question. Use of Equations (9) and (10) 
allows the construction of a true shear stress = rate of shear 


diagram from which the rheological class of the fluid can be determined. 
End Effects 


The equation for the rotational viscometer, Equation (9) 
is based on the assumption that the cylinders are of infinite length, 
i.e. that the viscous drag on the top and bottom of the inner cylinder 
is negligible compared with the drag on the cylindrical surfaces. This 
factor can be taken into account in one of the following ways: 
i) use of a very long inner cylinder of small radius with 
a small clearance between the inner and outer cylinders 
to make the end correction negligible, 

ii) calibration of the instrument with liquids of known 
viscosity, absorbing the end effect in the instrument 
constant, 

iii) reduction of the relative motion of the liquid above and 
below the ends of the inner cylinder by mechanical means, 
and, 

iv) measurement of the effect experimentally in order to 


include it in the calculations of the shearing stress. 
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Method i) imposes many difficulties in the design of a suit~ 
able apparatus. However, Shangraw et al (44) used this method with a 


ratio of length of inner cylinder to annular gap of 243 to l. 


Method ii) has been widely used and was used in this invest 
igation. However, it was found that the end effect varies with visco- 


sity. 


Method iii) includes the use of guard cylinders, which were 
introduced by Couette (4) and have been used by many investigators 
including Hatschek (23). Above and below the suspended cylinder, 
guard rings, (ise. cylinders of the same diameter), were rigidly 
fixed so as to be co-axial with it and to allow the minimum clearance 
between the free and fixed edges. Each of the guards is about 2/5 
of the length of the inner cylinder so that the viscous drag exerted 
by the base of the outer cylinder was not transmitted to the inner 
one. The liquid under examination filled up the interior of the suspen 
ded and guard cylinders. To prevent the liquid inside the inner 
cylinder from rotating, the opposite ends of the guard rings were 
filled with perforated plates. The required length of the guard cylin~ 
der varies with the annular gap and the velocity of rotation. In the 
apparatus of Couette, the guard cylinders had a length of 14 times the 


annular gap. 
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Method iv) has been widely used and is usually described in 
terms of an equivalent length of cylindrical bob. The plot of 
deflection versus bob length immersed, when extrapolated back is a 
straight line not passing through the origin, but intersecting the 
abscissa at a negative value equal to the effective end contribution. 
In further calculations, the measured depth of immersion may be 


adjusted by adding the equivalent length of the bob. 


Theoretically, end effect corrections of this type are 
suitable only for Newtonian fluids, since the consistency of a non= 
Newtonian fluid varies with the rate of shear and would vary across 
the base of the bob. However, in a well-designed apparatus, the end 
effect is small and the contribution to the total torque is small. 
Thus a relatively large error in the end effect correction itself is 


permissible. 


Stability of the Flow in the Viscometer Annulus 


In all the equations presented previously for the measure- 
ment of the rheological properties in a concentric cylinder viscometer 
it has been assumed that the flow is laminar. However, at high 
velocities, the flow may not be laminar and a suitable criterion must 
be used for the stability of the flow. A good treatment of the problem 
of stability of flow in the annulus of concentric cylinder viscometers 


is that of Chandrasekhar (2). 
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Probably the best known and most widely used correlation is 
that of G.I. Taylor who considered this problem from both the theoret- 
ical and experimental viewpoint. In 1923 Taylor (48) solved the 
equations of motion for the system of concentric cylinders containing 
a viscous fluid in the annulus, This allowed him to predict the 
onset of instability (turbulence) theoretically and his experimental 
work was an excellent confirmation of the theoretical work. Later 
Taylor (49) performed more experiments on different types of fluids 


and again the results were in excellent agreement with the theory. 


Taylor's theoretical work has been checked a number of times 
and found to be correct. Recent solutions of the theoretical problem 


include those of Meksyn (34), Chandrasekhar (3) and Di Prima (6). 


Experimental confirmation of Taylor's theory by other workers 
include the experiments of Lewis (29), Terada and Hattori (47) and 
Wendt (55). The most recent and probably the most accurate experiments 
are those of Donnelly (7) and Donnelly and Fultz (8). They check 
Taylor's theoretical predictions very closely. Taylor's working corr= 


elation is in the form of a log-log plot of Ut/y versus t/Ro where 


U = linear velocity of the outer cylinder 

t = gap between the cylinders 

y = kinematic viscosity of the fluid contained in the 
annulus 


Ro = -radius of the outer cylinder 
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From a knowledge of the geometry of the system it is possible to calcul- 
ate t/R, and obtain Ut/y from the correlation. From the value, the 


speed for instability to occur can be calculated. 


Another correlation, that of Schlichting (43) can be used 
to predict the onset of instability in a concentric cylinder viscometer. 


This correlation is a log-log plot of UR,/» versus t/Ro. 


In Appendix B, the results of calculations for the viscometer 
used in this work are given. It is seen that the agreement between 
the correlations of Schlichting and Taylor is very good. In this 
work, the velocity for instability in the annulus was always greater 


than 20 times the maximum speed of the instrument. 
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3. Literature Review 


A large amount of work has been done on the rheological 
behavior of suspensions of spheres. However, there is not complete 
agreement on the quantitative effect of the various physical para- 
meters. Good reviews of the rheology of suspensions have recently 
been published by Frisch and Simha (19) and Roscoe (41). In this 
section, a brief treatment of the most important studies is given with 
particular reference to the effect of particle size distribution 


upon the rheology of suspensions of spheres. 


Einstein (10), in a theoretical paper, solved the linearized 
Navier-Stokes equations for the case of a dilute suspension of spheres. 
The assumptions involved include the restrictions that the particles 
are rigid spheres their size being large compared to the size of the 
solvent molecules and yet small in comparison to the size of the 
apparatus in which the viscosity is measured. Another important rest- 
riction of Einstein's treatment is that the concentration of suspended 
spheres is infinitely dilute, so that the particles exert no forces on 
each other. With these assumptions, the expression obtained for the 


relative viscosity has the form: 


oes > uP = 1 + 2.5¢ (c ~» 0) C3) 
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Thus this case is a limiting result and states that the slope of the 


“We - ¢ curve has a value of 2.5 where it cuts the Y, axis. 


The Einstein equation was extended to higher concentrations 
in a theoretical work by Guth and Simha (21) who solved the Navier- 
Stokes equations by the method of successive reflections. The cal- 
culation was later revised slightly by Gold (22) and the result is: 


Me =toomuile ++) De Sel Lawes teat ee elt pee (14) 


Eirich and Goldschmid (12) found that Equation (14) 
held up to 6 volume percent for suspensions of glass spheres in 


mercuric nitrate - nitric acid solution. 


Other experiments by Eirich, Bunzl and Margaretha (11) 
were made with suspensions of glass spheres in mercuric nitrate - 
nitric acid solutions, mushroom spores in olive oil = tetrachloro= 
ethane mixtures, and yeast particles in water. The suspensions 
were Newtonian and the viscosity was measured in three types of 
viscometers - Couette, rising sphere and capillary viscometers. 

It was found that the Einstein equation was satisfactory up about 
5 volume percent and the following empirical equation represented 


the experimental results: 


3 
Ae = 1 fede t moc Paonia (15) 
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In 1948, Vand (50) (51) studied the problem of the viscosity 
of suspensions of spheres both theoretically and experimentally. In 
his calculation of the viscosity of a suspension, he allows for the 
effect of collisions between spheres as well as the hydrodynamic 
interaction. He obtains as a result: 


Wace 1 tatadeut Wease + Oh. (16) 


Experimentally, Vand measured the viscosity of suspensions 
of glass spheres in a solution of ZnI5 ~ water - glycerol and the 
results satisfy the empirical equation: 

2 3 
Vee te Deo ect glee Aen Oi ee) wi wa (17) 


in excellent agreement with his theoretical equation. 


Manley and Mason (31) found that there was an error in 
Vand's calculation of the collision frequency of the spheres and 


when corrected, Vand's theoretical equation becomes, 


Ne = | + 2,5¢ + 10.05¢> cae ae (18) 


Later experiments by Manley and Mason (32) using suspensions of glass 


spheres in ZnIy ~ glycerol gave as a result, 
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In a different approach to the problem, Mooney (35) considered 
the space crowding effect of the suspended spheres. This is the only 
theory to postulate an effect of size distribution upon the viscosity 
of the suspensions. For a suspension of groups of spheres, each group 


of different diameter, an equation of the form: 


m 
>) 
Inn, = 2.5 Drops (20) 
Asi 
i=l as jit 


is presented where A xi is the crowding factor. For a suspension 


containing only one size particle the equation reduces to : 


In = 2.5 ¢ (CAL) 
Ir (emis 


where k is {the self-crowding factor Ngae Mooney's approach is partly 
empirical in that the crowding factors must be experimentally determined, 


as they can only be predicted approximately. 
Effect of Size Distribution 


Ward and Whitmore (52) performed experiments with suspensions 
of spherical polymethylmethacrylate particles in lead nitrate - glycerol = 
water. The suspensions ysed were of varying size ratios (ratio of largest 
to smallest particle diameter) and the viscosities were measured in a 
rising sphere apparatus. It was found that the relative viscosity 
depended upon the size ratio of the suspension, decreasing with 


increasing size ratios 
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Eveson (13), using a Couette viscometer, studied similar 
suspensions to those of Whitmore and also used suspensions consisting 
of two widely different sizes. It was found that the unimodal suspen- 
sions were Newtonian fluids over the range investigated (0 - 20 volume 
percent), while the bimodal suspensions were pseudoplastic at concen= 
trations higher than 7.5 volume percent. No adequate explanation 
could be given for this behavior. However, it was found that the 
relative apparent viscosity was a minimum when about 50 percent by 


volume large particles composed the suspension. 


Another paper by Eveson (15) considered the influence of size 
distribution upon the viscosity of suspensions of spheres. Three 
size fractions of average diameter 378, 200 and 47 microns were used. 


The viscosity of suspensions of mixtures containing: 


50% 378 micron - 50% 47 micron 
25% 378 micron - 75% 47 micron 
75% 378 micron “ 25% 47 micron 
50% 200 micron = 50% 47 micron 
25% 200 micron ~ 75% 47 micron 


33 1/3% 378 micron ~ 33 1/3% 200 micron = 33 1/3% 47 micron 
was measured over a concentration range of 0 to 22.5 volume percent. 
It was found that for suspensions with a discontinuous size~distribution 
curve, the relative viscosity at a given solids concentration was never 
higher and was usually lower than the relative viscosity of the compon- 
ents of the mixture at the same concentration. This effect became pro~ 


nounced at volume concentrations greater than 17.5 volume percent. 
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Later, Eveson (16) investigated the effect of size ratio of 
the dispersed phase upon the relative viscosity of suspensions of 
spheres. There was disagreement on the effect of size ratio, Ward and 
Whitmore (52) having found a significant decrease in relative viscosity 
with increasing size ratio, while Eveson (13) and Higginbotham, Oliver 
and Ward (24) found little if any, effect of size ratio upon the 
relative viscosity. Eveson prepared eight samples, all having about 
the same average diameter, but with size ratios varying from 1.60 : l 
to 24 : i. He found that there was no effect of size ratio upon the 
relative viscosity. Eveson concluded that the effect of size ratio 
is unimportant and that a suspension possessing a discontinuous size 
distribution curve would be more likely to have a lower viscosity than 


a unimodal suspension. 


Some workers with synthetic rubber latices have studied the 
effect of size distribution upon the rheological properties of the latex. 
Synthetic rubber latex consists of small, almost spherical particles 
in a solution composed mainly of water. Fatty acids are used in the 
polymerization process and these are considered to be adsorbed on the 
surface of the latex particles. This adsorbed layer may be responsible 
for the pseudoplastic behavior of latices at higher concentration, this 
adsorbed layer being progressively sheared off as the shear rate 
increases, Maron and Madow (33) studied mixtures of two synthetic 
rubber latices, one with an average particle diameter of 980 R, the 


other with an average particle diameter of 1920 R. It was found 
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that the latices were Newtonian at concentrations below 25 volume per= 
cent and pseudoplastic at higher concentrations up to 55 volume percent. 
Also the consistency was a minimum at any volume concentration when the 


large latex was 75 = 80 volume percent of the total. 


Johnson and Kelsey (26) (27) carried out a practical study 
of the effect of size distribution upon the rheology of synthetic 
rubber latices. Since a Brookfield viscometer was used and the 
results for the pseudoplastic latices interpreted as if they were 
Newtonian fluids, their results can give only a qualitative picture 
of the effect of size distribution. ‘They used three latices of average 
particle diameter 950, 1710 and 3250 £ and found that when the large 
particles comprise about 75 volume percent of the total and when the 
difference between the small and large diameters is large, the 


greatest reduction in "viscosity" was noted. 


Sweeny and Geekler (45) studied the rheology of suspensions 
of glass spheres in ZnI5 - glycerol ~« water. They found the suspen- 
sions to be pseudoplastic and characterized the rheological behavior 
by the terms viscosity at low rate of shear !) y and viscosity at 
infinite rate of shear co As part of this study they investigated 
the effect of size distribution on the rheology of the following 
suspensions at 55 volume percent. 

i) 100% 262 microns 


Li) 754 262 microns - 25h 164 microns 
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25% 97 microns 


iii) 75% 262 microns 
iv) 75% 262 microns - 25% 36 microns 
v) 75% 262 microns - 25/4 13 >microns 
It was found that the relative viscosity at zero rate of shear and at 
infinite rate of shear decreased from suspension i) to suspension v) 
in the above table indicating the influence of the diameter ratio 


upon the consistency of bimodal suspensions. 


In an unpublished paper, Sweeny (46) recently extended the 
previous work to include results of suspensions of glass spheres. 
Measurements were made at 40 and 55 volume percent only and four 
diameter ratios were considered as in the previous paper. However, 
measurements were made at 25, 50 and 75% of the larger spheres at 
each of the four diameter ratios. It was found that the minimum 
"viscosity'' occurred when about 70 volume percent of the total 
solids consisted of large spheres. The diameter ratio also had a 
significant effect upon the observed viscosity, the maximum reduction 
in viscosity occurring with mixtures of the largest and smallest 
spheres. Mooney's equation (Equation (20)) was found to represent the 
data well, but there is an error in the calculations which invali- 


dates this conclusion. 


From the published literature the following conclusions 


can be drawn: 
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i) There is no general agreement among the equations of 


various workers relating the dependence of viscosity upon concentration. 


ii) Suspensions of true spheres in polar media usually behave 
as Newtonian fluids. Thixotropic and pseudoplastic behavior has been 
noted. In some gases, these anomalies can be explained in terms of 
surface effects (adsorption), or flocculation of the spheres. Within 
the experimental error, no dilatancy has been found in suspensions of 


spheres. 


iii) The particle diameter of the spheres in a unimodal 
suspension has no effect upon the viscosity over a wide range of 


diameters (10 to 1000 microns). 


iv) When the solid phase has a continuous distribution of 
diameters, there is little effect of size upon the rheological behav- 
ior of the suspension over a large range of size ratios (the ratio 
of the size of the largest sphere in the suspension to that of the 


smallest) « 


v) When the solid phase has a discontinuous size distri- 
bution, the consistency of the suspension is lower than that of the 
unimodal components at the same volume concentration, This effect is 
pronounced at volume concentrations exceeding 20 volume percent and 


is greatly influenced by the diameter ratio of the two components in 
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a bimodal suspension, the suspension having the highest diameter 


ratio also having the lowest consistency. 


vi) The effect of size distribution has b@en carefully 
studied up to concentrations of 20 volume percent. Above this data 
are available only at 40 and 55 volume percent. The need is evident 
for work at intermediate concentrations to tie together the effect of 


size distribution upon the rheology of suspensions of spheres. 
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ITLL, EXPERIMENTAL EQUIPMENT AND METHOD 


1. Viscometer 


A concentric cylinder viscometer was used to determine the 
rheological properties of the suspensions in this investigation. A 
"Fann V=G Meter'' (Model 35) manufactured by the Fann Instrument 
Corporation, Houston, Texas previously used by Lindley (30) .and 
Aleorn (1) in their work on the rheology of suspensions was also 
used in this work. Figure 3 shows the basic features of the instru- 
ment. When the sample cup is raised, fluid enters the annular space 
between the two co-axial cylinders. The rotating cylinder or cup is 
driven through a system of gears by a 115 volt = 50 cycle dual 
speed synchronous motor. The torque exerted by the viscous drag on the 
inner cylinder or bob is balanced by a double~wound helical torsion 
spring. One end of the spring is coupled to the upper end of the 
bob shaft; the other end is secured by a mounting clamp. The angular 
deflection of the bob is read from the scale on the upper surface of the 
indicating dial which is mounted on the bob shaft. The scale is 
marked at l-degree intervals through 270 degrees and is numbered 
at 10-degree intervals. An optical reticule is provided for accurate 


reading of the scale. 
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FIGURE 3 


PHOTOGRAPH OF THE VISCOMETER 


In the standard model, cup speeds of 3, 6, 100, 200, 300 
and 600 rpm are obtained smoothly and rapidly without stopping rotation 
by means of a three position gear shift in conjunction with the two 


position switch controlling the motor speed. 


The instrument used in this work was modified in three ways. 
Firstly, the stainless steel cup used on the standard model was replaced 
by a copper cup with an annular jacket. Ethylene glycol solution from 
a large (10 Imp. Gal.) constant temperature bath was circulated through 
the sample cup jacket by means of an immersion pymp. The temperature 
in the sample cup was held at 20 + ule by means of a Thermotrol 
(Hallikainen Instruments, Berkeley, California) controlling the temper- 
ature of the bath. Secondly, additional shear rates were obtained by 
the interposition of a 9 : 1 reduction gear between the synchronous 
motor and the body of the instrument. This modification permitted 
greater coverage in the region from 6 to 100 rpm. Thirdly, two 
additional cups of larger diameter than the standard cup were machined 
from brass. This permitted the use of large particles in the suspensions 


without binding or wall effects. 


In addition to the standard helical torsion spring, two 
weaker springs were purchased from the manufacturer. These springs were 
1/2 and 1/5 as strong as the standard spring, and allowed measurements 


to be made on fluids of low viscosity. 


el ee oe 


O0¢ ,0OS ,OOL .d ,£ 26. abseqe quo .{sbom busbaate edt al 
moizasox gaiqqese jvodtiw yibiqat bas <idsoome bentstdo sxe aq 002 bas 
owt edt diiw notiomstmos at a¢ide 29g moltinoq serntt a 19 canen xd 
«beege tosom si3 galiloximes dotlwe noltieeq 


-oyow souls sit GURRNDG sme few ald at a 
bsosiqes saw Iobom busbnate edt no besw gro lenge exeint 1 os eeltextt 
moxs soltulan Looyly emetyds’ stodnnt alam A ib son ct 
dgvowdt bstsluouks esw diad oxstexsqms tamganoa (180 OL) sgxal 
omieteqmes edt -;494 molesamni =» 20 amnom ys aednat quo signe ads > 
JortoumaaT a 2 aman yd 9°10 + OS 20 Breit aa que afman oft at 
-xeqmad edi galiloxtuos (atarohtiad cyofedaol ysamemnseud anmtnststalty 
<i benteddo otew aetex uaeds [asotiibhs «yibaooe® vdted alt 3o ome tan 
sven oe inet Se eset ae : a 
hesa.iermeq sottsokiibom «tlt end | 
ows ,vlbuldT ane OOL oF 0 aoxt 
bentrioan oxmw quo bimbmate est mtd | . 
cootacngeve ads seins sale So can fot Milian aie ssenrd moet 
-staatte aw vo wo gathadd suadity 
a. 
ows ,yituqe molowed Meotied bmabante eds of molsibba el 


eTaw agabxqe seedT -sevetashamen at cot heetheiey/seew Sgiine| Seine 


’ 


aJooon resem bewolls bas ygataye SUNN 5.08% Sa : va 


out vr _ 


*, | = ae 


a 
| 


2 


npn he = 74 eee 


The various cup~bob combinations were calibrated with 
Newtonian standard oils. The viscosities of these oils had been 
determined at 70°F in Master Cannon~Fenske capillary viscometers. 
Four oils were used, of viscosity 19.1, 56.7, 93.6 and 147.9 


centipoise. The results of these calibrations are presented in 


Appendix A. 
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2. Suspensions 


The suspensions used in this work consisted of spherical 
particles of styrene~divinylbenzene resin in an equal density solution 
of glycerol, isopropanol and water. The styrene-divinylbenzene 
copolymer resin (8% divinylbenzene) was obtained from Dow Chemical 
Company, Midland, Michigan and three nominal sizes were used. They were 
40/60, 100/200 and 200/400 mesh (U.S. Standard Sieve Series). 
Duplicate sieve analyses employing both dry sieving and wet sieving 
with isopropanol were performed using stainless steel U.S. Standard . 
Sieves. ‘The results of these analyses are given in Appendix D and 
they are plotted in Figure 4 to show the size distribution of the 
original materials. The density of the beads was determined by the 
standard -pycnometer method and the average density was found ta be 


1.0487 gm/ml. 


Samples of three powders were photographed through a 
microscope (71X) to give an indication of the size distribution and 
shape of the particles. These photomicrographs are shown in 
Figures 5 and 6. It can be seen that the particles are all spheres and 


that the size distribution for any powder is not large. 


Suspending solution No. 1 had the following approximate 


composition: 
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FIGURE 5 


PHOTOMICROGRAPHS OF THE 40/60 MESH AND THE 100/200 MESH POWDERS (71x) 


== 


FIGURE 6 


PHOTOMICROGRAPH OF THE 200/400 MESH POWDER (71X) 
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Weight Percent 


Glycerol 40 
Lsopropanol yaa) 
Distilled water 35 

100 


The density and viscosity of the solution were determined at 20.0 C, 
The standard pycnometer method was used in the density determination 


and the viscosity was measured in a Master Cannon-Fenske capillary 


viscometer. 
Density ~ 120491 gm/ml at 20.0 C 
Viscosity - 12.14 centipoise at 20.0°C 


Near the completion of this investigation suspending solution 
Noe 1 was exhausted and a second solution was prepared. It had the 
following composition and properties: 


Weight Percent 


Glycerol 45 
Tsopropanol 30 
Distilled water 25 

100 
Density - 1.0494 gm/ml at 20.0°C 


fe) 
17.29 centipoise at 20.0 C 


Viscosity 


eh2< 


3. Auxiliary Equipment 


Temperatures were measured by a sixteen inch mercury-in- 
glass thermometer (Scale: - 1° to 51° C). The thermometer was 
calibrated against a standard platinum resistance thermometer. The 
calibration data and curve are given in Appendix C. All temperature 
readings were made with a thermometer reading lens which allowed 


the thermometer to be read to 0,.02°C. 


Other equipment used included an analytical balance accurate 
to 0.001 grams with a maximum load of 500 grams and a graduated 
cylinder marked at 1 milliliter intervals. These two pieces of 
equipment were used to determine the volume concentration of the 


suspensions. 
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4. Experimental Method 


Two procedures were followed in measuring the rheological 
properties of the suspensions. In the first, a known volume of 
suspending solution was poured into the sample cup and allowed to 
attain the bath temperature. When the temperature, measured with the 
standard thermometer, was 20.0°C + 0.1°C, the sample cup was raised 
so that the fluid covered the top of the tape on the viscometer cup. 
The viscometer was switched on the lowest speed at which a reading 
could be obtained. After the stable reading had been noted, the vis- 
cometer was switched to the next highest speed. This was continued 
until the highest speed was reached and the reverse procedure was fol- 


lowed, obtaining the readings at successively lower speeds. 


Next, solids from a weighed beaker were poured into the 
sample cup and the beaker reweighed. The sample cup was placed on 
top of a magnetic stirrer and stirred by means of a Teflon~covered 
magnet for about 2 minutes. The temperature of the suspension was 
read and the dial deflections versus speed obtained as before. In 
order to avoid cumulative errors in weighing, the spheres were 
added from a full beaker and the beaker weighed when full and after 


each addition of solids. 


This procedure was followed until the torque was too great 


to be measured with the weakest spring (No. 3). Then the No. 2 spring 
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was inserted and the readings continued. Eventually, above 55 volume 
percent, the suspension became so viscous that readings with the 
strongest spring could be obtained only at very slow speeds and the 


readings became so erratic that at this point the test was terminated. 


Although this procedure would seem to be satisfactory, a 
problem arose in the calculation of the volume concentration. Over 
the course of a run ( 4 to 6 hours) the viscometer cup was out of the 
sample cup for long periods of time, Evaporation of solution occurred 
and panridants of the suspensions were lost when the thermometer was 
removed after a temperature reading. Also since the sample cup was of 
fixed volume and enough solution was added initially to allow a reading 
to be obtained, after some solid had been added, the sample cup was 
full. In the process of removing suspension between tests errors may 
have occurred in the volume withdrawn. These problems were brought out 
by samples taken during a rune It was found that the measured concen 
tration was higher than that calculated from the initial volume of 
solution plus the solids added minus the amount of suspension removed. 
This was also checked by experiments made at an accurate volume concen= 
tration, by mixing together the amount of solid and liquid necessary 
for a particular desired concentration. On the basis of these two 
cheeks, the calculated concentrations were corrected at the higher 
volume fractions; however, the correction was not a constant and 
increased as the concentration increased, For this reason a second 


procedure was adopted for later runs. 
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In this method, an amount of solution required to produce 
a suspension of the desired solids concentration and a total volume 
of 240-250 ml was poured into the sample cup. The amount of solids 
necessary to give the desired solids concentration was weighed out 
by means of the analytical balance and the suspension mixed together 
in the sample cup with the magnetic stirrer. The temperature of the 
suspension was measured with the standard thermometer; then the sample 
cup was raised and the dial deflections read as the speed was increased 
to maximum and decreased to zero. The sample cup was allowed to 
remain around the viscometer cup and bob for 5 minutes. Then the 
dial deflection~speed sequence was run as before. The sample cup 
was lowered, the sample discarded and the viscometer cup and bob and 
the sample cup cleaned with water in preparation for a new test. 
In this way, the solids concentration was known accurately and the 
major portion of the error is in the viscometric measurement. 
However, this method requires much more solid and solution than the first 


procedure. 
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5e Mechanical Difficulties 


Bearing Friction 


In preliminary calibration experiments, it was found that 
dial deflection versus speed data for Newtonian standard oils did not 
fall on a straight line with a slope of wnity when plotted on log-log 
paper. There was a definite curvature at low speeds, the data being ' 
higher than the line with a slope of unity. It was found that this 
was a result of friction in the bearings separating the bob shaft 
from the rotating sleeve that drives the cup. It was suggested that 
the bearings were old or worn; however, replacement with new bearings 
gave the following results. If the viscometer dial was set to zero 
when the cup was not in motion, a deflection of 1 to 2 degrees was 
observed when the viscometer cup was removed and the instrument 
turned one This deflection accounts for the higher readings at low 
values of speed. For example, at low speeds with a deflection of say 
10 degrees, 1 to 2 degrees is an appreciable fraction of the total, 
while at high speeds with a aes ection of say 220 degrees, this amount 
is almost insignificant. It was found that the error due to bearing 
friction could be partially eliminated by setting the dial deflection 
to zero while the viscometer was running in air without the cup 
attached. This was done at a speed of 300 rpm and it was found that 
at very low speeds the dial deflection at no load was slightly 


negative, while at 600 rpm, the deflection was about 0.5 degrees. 
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This fact limited the precision of the dial readings for the instrument 


to + 0.5 degrees. 


Temperature Variations 


Although the bath temperature control was wery good, problems 
were encountered whenever the ambient temperature differed from the 
bath temperature by more than about 0.5°C. Before a test, the viscometer 
is at the room temperature. When the sample cup is raised, covering 
the viscometer cup and bob, any difference between the ambient tempera~ 
ture and the sample cup temperature results in a temperature gradient 
along the viscometer annulus. The viscosity of the suspending solution 
changed about 6% per degree centigrade change in temperature, so only 
a slight difference between the room temperature and the bath temper= 
ature had a significant effect upon the measured rheological properties. 
For this reason, in all tests, the room temperature was adjusted to 
7050.-10,5226 eens any tests were performed. It would seem that 
a constant eipevarune room or cabinet is necessary for ease in 


taking measurements. 


Another problem in temperature control is the heat buildup 
upon prolonged shearing. It was found that after shearing a Newtonian 
standard oil for 2 to 3 minutes, at 300 rpm, a slow decrease in dial 
deflection with time was noted. This was probably due to inefficient 


conduction of the heat produced in shearing. In this viscometer heat 
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produced in the annulus must first be conducted through the viscometer 
cup to the bulk sample of suspension. Thus another temperature gradient 
exists, after shear, between the wall of the sample cup and the annulus 
of the viscometer. This inefficiency in removal of heat produced 

in shearing meant that accurate viscometric measurements could not 


be made continuously for more than about 2 minutes. 


Cylinder Alignment 


It was noticed that the viscometer cups or outer cylinders 
did not turn concentric with the viscometer bob, but precessed as 
much as 3.5 thousandths of an inch. It is to be expected that this 
misalignment would increase the dial deflection since the minimum 
dissipation of energy would occur when the cup and bob are concentric. 
It was considered, however, that errors due to this cause were 
insignificant compared with those due to temperature control and 


other mechanical difficulties. 
Oscillation of Viscometer Bob 


In preliminary tests to determine the suitability of 
suspending solutions, tests were made with solutions of viscosity 
3.2 and 6.5 centipoise. However, it was found that the dial readings 
oscillated widely with the 6.5 centipoise solution for about 1 minute. 


With the 3.2 centipoise solution no stable reading was obtained. For 
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this reason, the 12.14 centipoise solution used in this work represents 
a practical minimum with which accurate viscometric measurements are 
possible. Weltmann (54) has stated that in order to obtain dial 
deflections rapidly, the system should be slightly overdamped. In 

the case of the Fann viscometer, the heavy metal bob ensured that the 
system is underdamped unless fluids of viscosity 50 centipoise or more 
are studied. For accurate work on fluids, the viscosity of which is 


as low as water, a light plastic bob with low inertia would be required. 


Sample Level in Viscometer Annulus 


In preliminary calibration runs with a Newtonian standard 
oil, it was found that dial readings at 600 rpm were consistently 
lower than those required for a straight line of dial deflection 
versus speed. An insight into this phenomenon was obtained when it was 
found that if the sample cup were not raised to the scribe mark on 
the viscometer cup, a linear response was obtained including the 
600 rpm speed. The article by Savins and Roper (42) explained this 
matter fully. These workers used a Fann VG Meter and found the deflec- 
tion at 600 rpm to be 1.5% low compared to the straight line drawn 
through zero and the 300 rpm reading. Whereupon, they constructed a 
transparent plastic viscometer cup and a transparent sample cup in 
order to observe the flow pattern in the annulus. When a small amount 
of dye was introduced near the base of the bob, it was observed to flow 


upwards in the annulus and out the ports in the upper portion of the 
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viscometer cup. Also the position of the liquid level in the annulus 
was found to depend upon the cup speed. These workers found that, 
" at the 300 rpm setting the liquid level dropped appreciably 
below the rest level but very little of the upper portion of 
the bob was exposed above the liquid. The drop in level at 
the 600 rpm setting was so pronounced that in most cases the 
conical upper surface of the bob was completely exposed 
above the sample. However, by sealing the sleeve ports 
with tape, the circulation was eliminated with the result 
that the liquid level did not drop below the upper tip 
of the conical portion of the bob at either rotor speed." 
Since this was reponsible for the apparent low readings at 600 rpm, 
the sleeve ports were sealed with electrical tape for all subsequent 
calibration and experimental work. With this technique the 600 rpm data 
lay on a straight line through the origin and the lower speed data 


points. 


Apparent Thixotropy in Sphere Suspensions 


In preliminary experiments with polystyrene spheres (Dominion 
Anilines and Chemicals, Toronto), thixotropic behavior was noted at 
concentrations as low as 20 volume percent. That this was thixotropy 
and not a temperature effect was evident from the rapidity with which 
the dial deflection decreased in the first minute of shear. After this 
period the dial deflection decreased slowly making it difficult to sep- 


arate the effects of temperature increase and thixotropy. 


avivhes ¢d3 mt fever sant as weg i bet A 


hi 


tad? bevel ssedae < smite toons qe mesloe wing 
= . yidats saqgs basgoxs Lowe bispt! ‘ety guts 


‘a 
a : to moltaog se0n oda to sizat yw , tind | ve eon oid 
8 +e "i 


jn fovel al qoexb eit ehiapit aits svete Banos em So a 
/ sit gsefe J2som at gade beomeowery of ast gnisise -_ “008 oar 
f : >. : 6 - Hn re -* 
| bezoges yledelques ssw dod srt to eonttéa seqay [eeknoo 


re | ae 2 SS 
atsoc Sveeia ena gnilacs yd ,9veWwoH ooh ORE ovyods 
i : 7 _ aul 


sivast eft Atiw bejsntalle caw coltalnonts ome ogad 3 tw 


7” 


: mat 
qks zeqqu ed? woled goxb som bib Loved biophi ond anda 
| & 
" baeqa rogox stedtle dod ad? bs fy ot3 20g [notes sdt x0 
ra s 4 7 : 
4 ogy 900 Js eagnibnex wel trexzhoqus ef tot sidisgoqs: saw =e e 


Ireupsadye iis tot eqka San) wtoals “tiw balsea oruw 81x04 ovesta 

siah mgqz 000 em? suphadoe? atdy W . strom [st coed tegxs bos wot wal | 
siab hesqe tewol add bre atgizo oid quo wnt sighessa & 

bir ruil sezaac! 
| (Men awa 
T a 
7 nohileod) se redge esetyeeylon atiw at: wn) egies _ in 
S ? je Sethe avy volver stowsd oxiay ‘ ane ata + Lawl 

a vy yqortoxldt nsw eht taaT tiene omulov O° r woh Cw ei st 


ry 


a ’. = 


| alow dolw yaibieas oity ih CL. eae — toette | a mt 
i. 7 7 : “= : 
_ ali ‘wen dA « THe (16 th etuale vari ety me 5 oh rz 

;, —, - | fn 
“geo Of Jloollith 2) gabe eieele ~ 


7 7 a 
stqowonldd | pom d nn 4 dais 


5 
a 
‘Was Lelia ae 


Since it was not desired to study the thixotropic breakdown 
of a suspension and also because it is difficult to envision a suspen- 
sion of spheres as a thixotropic fluid, attempts were made to disperse 
the spheres with chemical dispersants. However, a large number of 


dispersants of various classes of compounds failed to have any effect. 


At this time, the styrene-divinylbenzene copolymer spheres 
were received and in experiments with these particles, thixotropic 
behavior was noted at concentrations greater than about 30 volume perc~ 
ent. Thus both types of spheres appeared to be similar in their sus- 


pension properties. 


The reason for the thixotropic behavior is thought to be 
flocculation of the spheres. In the process of shear, these 
ageregates are broken down, decreasing the consistency with time as 
shearing progresses. Since flocculation was considered to be the cause 
of the problem, the forces involved which create the aggregates were 


considered. 


Three possible types of forces can cause flocculation. These 
are chemical forces, electrical forces and physical forces. Firstly, 
the thixotropy did not seem to be a result of chemical forces, in 
that dispersants of various chemical types had no effect. Electrical 
(electrostatic attraction) forces also did not seem to be the answer 


because the solution is polar and charge induced on the particles 
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should slowly drain off, especially since the suspension was placed 
in a copper cup which was grounded. Another reason why electrical 
charges were not considered to be the cause of the thixotropy is that 
addition of ionic compounds such as NaCl and. Pb(NO 3) 9 had no effect. 
The addition of compounds such as these increases the conductivity of 
the solution and would tend to dissipate electrical charges very 
quickly. Physical forces, in particular surface tension, appear to 
be the most plausible explanation of the thixotropic behavior. The 
particles are originally in air and when placed in the suapenciing 
solution, they would presumably retain some air on the surface. This 
air bubble could hold two or three particles together by the suface 
tension between the solution and the air bubble. Since these forces 
are weak, shearing breaks up these doublets and triplets. However, 
the air bubble would not necessarily rise as it is very small and 

the head of solution over it may effectively prevent its escape. 

Some evidence in support of this explanation may be found in the 
determination of the density of the plastic particles, In this 
determination, the particles were placed in a density bottle and a 
heavy oil poured over them. Vacuum was applied to the bottle to 
remove all traces of air. It was found that small bubbles of air 
continued to rise from the solid for as long as four hours on the 
vacuum pump. From this it can be seen that small traces of air 

might never come to the surface of the suspension. Unfortunately, 
attempts to apply a vacuum to the suspension caused the solution to 


boil, because of the large amount of water and isopropanol which it 
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contained. The experience of Vand (51) also supports this explanation. 
He used a Couette viscometer and stirred the suspension of glass spheres 
by bubbling air through it. It was found that atconcentrations above 

40 volume percent the dial deflection immediately after air stirring was 
much greater than normal, and decreased with time exponentially to its 


normal values 


Since the true effect of the spheres in the liquid is noted 
when the aggregates are completely broken, mechanical shear outside 
the viscometer was attempted. The sample cup was placed on a magnetic 
stirrer and a Teflon-covered magnet placed inside the cup. It was 
found that shearing by this method gave no further change in the rheol= 
ogical properties after 2 minutes and the suspension remained in a 
broken down state for up to 30 minutes after shearing ceased. The 
procedure adopted for viscometric measurements involved 2 minutes of 
shear with the magnetic stirrer immediately before a test was made. 
Tests were conducted on a particular suspension for up to 20 minutes 
after shearing by this method. By this standard procedure, a sample 


upon which reproducible measurements could be made was prepared. 
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EXPERIMENTAL RESULTS AND DISCUSSION 


Viscometric measurements were made on the following suspensions 


over a concentration range of 0 to 50 volume percent: 


i) 
it9) 
iii) 
iv) 
v) 
vi) 
vii) 
viii) 
ix) 
x) 
xi) 


40/60 mesh 
100/200 mesh 


200/400 mesh 


- Material A 


- Material B 


- Material C 


50% A - 50% B 
254A = 75% B 
75%ZA = 25% 8B 
50%ZA = 50% C 
D5 k= BTS 
754 A= “95% € 
50% B = 50% C 
Bo 1 Se A He SBF /a% B= 33 -1/37"¢ 


It was found that the suspensions behaved as Newtonian fluids over the 


entire range 


of concentration studied. 


The basic experimental data 


are given in Appendix E, and the calculated results are given in 


Appendix F. 
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1. Unimodal Suspensions 


The relative viscosity~concentration data for the three 
unimodal suspensions are plotted in Figure 7. It can be seen that 
there is a significant difference between suspensions composed of 
material B as received and those composed of materials A or C. This 
difference did not appear to be a dependence of viscosity upon particle 
size as material B is intermediate in size. Nor did it appear to be a 
result of size distribution, since all three original powders had 
essentially the same narrow distribution of particle sizes. Only one 
physical difference between the powders was noted. Samples A and C 
were free-flowing, that is, the angle of repose in air was very small. 
On the contrary, sample B was not free-flowing. Material on a spatula 
would not fall off even when the spatula was held vertically. The 
reason for this difference in dry flow properties was Luought to be due 
to solvent which caused adherence of the particles. This solvent was 
presumed to be the material in which the particles were polymerized, 
and was not completely removed when the particles were dried. For this 
reason the particles were cue isopropanol to remove the solvent. 
About 2 pounds of material B was mixed with isopropanol to form a dilute 
slurry. A portion of this slurry was poured on a 400 mesh sieve and 
the isopropanol allowed to drain off. Fresh isopropanol was then 
poured over the sieve to remove any trace of solvent remaining. 
The sieve was set to dry in an oven at 40°C. When dry the powder 
was found to be free-flowing. This procedure was continued until all 


the slurry had been washed in isopropanol and dried. Figure 7 shows 
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that the washed material exhibits a viscosity~concentration relationship 


similar to materials A and C. 


An interesting observation is the range of concentration over 
which the Einstein equation is valid. A plot of the observed relative 
viscosity divided by that calculated by means of the Einstein equation, 
Equation (13), is shown in Figure 8 for unimodal suspensions. It can 
be seen that the Einstein equation would be correct to within 1% only 


at volume concentrations up to 2 volume percent. 


An attempt was made to fit the data to a power series in 


concentration of the form 


3 4 


Mx = lt 2Sce+actt bev tect . . (22) 
by the method of least squares, The coefficient 2.5 was chosen for 
the second term because theoretical and experimental work has shown 
this constant to be correct. Also the curvature of a graph of ae = 
concentration is such that accurate estimates of this constant can 
only be made from measurements at volume concentrations of less than 
5% and this requires very aceurate determination of the relative 
viscosity. Since only a few data were obtained in this concentration 


range, the value 2.5 was taken for the first term of a power series in 


concentration. 
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Since the relative viscosity increases continuously and mono- 
tonically with concentration, at -/de? > 0, and all coefficients in 
Equation (22) must be positive. However, using the method of least 
squares to fit the data, it was found that negative coefficients were 
obtained in some cases. Also, coefficients that were negative when 
say a third power series was tried, became positive when the series 
was changed to a fourth power. For these reasons it was concluded 


that fitting by the method of least squares was not satisfactory. 


A second method was tried. This consisted of plotting 
(N, - 1)/e versus c to obtain the coefficient of ¢ as the inter- 


2 as the slope of the linear 


cept at c = O and the coefficient of c 
portion of the graph. This plot for unimodal suspensions is shown in 


Figure 9 and it can be seen that a power series of the form 


ls = ee de 10 ee (23) 
is indicated. It was found that the higher terms of this series could 
not be determined accurately by this method. The range of fit of 
Equation (23) is shown in Figure 10 in which the observed viscosity 
divided by that calculated by means of Equation (23) is plotted as a 
function of concentration. It can be saen that the equation is correct 
to within 5% at concentrationsup to 14 volume percent. Equation (23) 
checks the theoretical equation of Manley and Mason, Equation (18), 


very closely. 
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It can be seen that there is little difference among the 
results for unimodal and the various bimodal suspensions at concen- 
trations below 20 volume percent. Thus all the results may be fitted 
by Equation (23) at low concentrations. For this reason, a method of 
correlation was sought which would show the effect of size distribution 
at high concentrations. The method of Vand (51) was used. In this 
method, the relative fluidity, 0,5 is fitted to a power series in 


concentration of the form 


$.-= 1 - 2.5¢ + act + be? (24) 
Because the curve of 9. - c is much more nearly linear than the curve 
of ae ~ c, only three terms need pe used. If the equation is recast 


in the form 


Wes 25 = -~a -—be (25) 


a plot of the left hand side versus c should give a straight line of 
slope -=b and intercept «a at c = 0. It was found that this method 


worked well and was used for all of the other data on the bimodal 


and trimodal suspensions. The calculations and plots used to obtain the 


coefficients a and b are given in Appendix H. 


For unimodal suspensions, the fitted equation is 
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In Figure 11, a test of the fitted equation is given by plotting 

“x (observed) X 6. (calculated from Equation (26)) versus concentration. 
It can be seen that most of the points are within + 5% of unity. This 
method of correlation is satisfactory for 0, because of the relative 
linearity of the 6. - c plot; the series can then be truncated at 

3 terms. However, if this three term series is converted to a series 

of “re - c by long division, the resulting series is only very 

slowly convergent and cannot be truncated. For this reason, all of 


the fitted equations were left in the 6, form. 
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2« Bimodal Suspensions 


The data for bimodal suspensions of A and B are given in 
Figure 12. The line for the unimodal suspensions is given for reference. 
It can be seen that there is little difference between the results for 
50% A = 50% B and those for 754A = 25% B. The data were fitted 
by the same procedure used for unimodal suspensions and the resulting 


equations are: 


25h A. = 715k B 


g. = 1- 2.5¢+ 0.2c7 + 1,96c2 (27) 
50%.A = 50%.B and 754A » 95% B 
$6. = 1-2.5¢+ 0,.8c~ + 0.922 (28) 


These equations were tested by the same procedure used for unimodal 
suspensions and the results are given in Figures 13 and 14, Most of 


the points are within + 3% of unity in Figure 13 and + 5% of unity 


in Figure 14. 


Data for bimodal suspensions of A and C are presented in 
Figure 15. A comparison with Figure 12 shows that the viscosity at 
any concentration is slightly lower than that for bimodal mixtures of 


A and B. 
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The equations pbtained by fitting the data to a power series 


in concentration are: 


254 A - 75% C 


$. = 1 = 2.5¢ + 0410c” + 2.26c° (29) 
50% A = 50% C and 75%ZA = 25%6 
G. = 1 2.5c + 0.97 + 0.915c° (30) 


Figures 16 and 17 show the results of the test of these equations against 
the experimental data. Most of the data points are within + 5% of 


unity. 


Data for the bimodal suspension of 50% B and 50% G and for 
the trimodal suspension of 33 1/3%4 A = 33 1/34 B = 33 1/3% C are 
given in Figure 18. Since the data for suspension of 50%Z B = 50Z C 
differed only slightly from that for unimodal suspensions, the other 
combinations, 25% B - 75% C and 75% B = 25% C were not investigated. 
It may be said that the three component suspension is not really 
trimodal since the diameters of B and C are very close. Thus this 
suspension would be expected to behave as 33 1/3% large ~ 66 2/3% 
small particle suspension. The fitted equations for these data are 
as follows: 


f. = 1 2y5e t Oc* + 2,4e3 (31) 
33. 1/97, Ao 33 1 /G Box 33. 1/32 6 
$B. = 1 = 2.5¢ + 0.5c7 + 1.6e° (32) 
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Figures 19 and 20 are plotted to test the fit of these expressions 
to the experimental data. The data points are + 3% of unity in 


Figure 19 and + 5% of unity in Figure 20. 


In Figures 21 and 22 respectively, the coefficients of 
c* and c3 in the power series in G.. are plotted as a function of 
diameter ratio of the larger to the smaller spheres in the suspension.« 
The diameter ratio is the volume-surface mean diameter of the larger 
particles divided by the volume-surface mean diameter of the smaller 
particles. This diameter was chosen to make the results consistent 
with those of Sweeny (46). The point for the 33 1/3% larger particles 
is the trimodal suspension, the value of d/D obtained by considering 
the average of the diameters of Materials B and C as one size and the 
diameter of Material A as the other. It can be seen that there 
isa fair correlation between the coefficients of the power series 
and the diameter ratio. The coefficient of o2 appears to increase 
and the coefficient of ec? to decrease with increasing diameter 


ratio and seem to approach a constant value at high diameter ratios. 


The data for the bimodal suspensions have been plotted as 
the ratio of the bimodal relative viscosity to the unimodal relative — 
viscosity under the same conditions. The ratio is plotted as a 
function of percent larger particles. with volume percent total solids 
as a parameter in Figures 23,24 and 25. Figure 23 contains the 


results for bimodal suspensions of B and C, for which the diameter 
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ratio, the ratio of the average diameter of B to that of C, is 1.8. 
Similarly, Figure 24 presents the results of bimodal suspensions of 
A and B, for which the diameter ratio is 4.5 and Figure 25 presents 
the results of bimodal suspensions of A and C, for which the diameter 
ratio is 8.25. It can be seen that there is a pronounced effect of 
diameter ratio, the ratio of the relative viscosities for the large 
diameter ratio being much lower than that for the small diameter 


ratio. 


Figure 26 compares the present work with the results of 
Sweeny (46) at 40 volume percent, the only value at which Sweeny has 
results which are also available from the experimental data. The figure 
is a plot of the ratio of the relative viscosities as a function of 
the volume percent larger particles with diameter ratio d/D as a 
parameter. In general, the agreement between Sweeny's results and the 
present work is very good. Sweeny's results for a d/D of 7.3, 
however, lie in the same region as experimental data from this work 
Eon a d/D Of S25. tm fehis work all references to Sweeny's data are to 


viscosities at zero rate of shear. 


Cross plots of Figures 23, 24 and 25 were made at 25%, 50% 
and 75% of the larger particles in a bimodal mixture. The ratio of 
the relative viscosities is plotted as a function of diameter ratio 
with total volume percent solids as a parameter in Figures 2/7, 28 
and 29. On these figures are plotted data of other workers including 


Eveson (15), Sweeny (46) and Sweeny and Geckler (45). It can be seen 
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that the agreement between the data of other workers and the present 
work is generally good, although the results of Sweeny at 55 volume 
percent solids and with bimodal suspensions containing 50% of the 

larger particles lie on the same curve as data for 50 volume percent 
solids obtained in this work, It can be seen that there is a pronounced 
effect of diameter ratio, especially at higher concentrations of 

solids, the higher the diameter ratio, the lower the relative viscosity 


compared to unimodal suspensions. 


An explanation of the effect of diameter ratio upon the 
ratio of the relative viscosity of bimodal suspensions to the 
unimodal suspensions may be found in consideration of the packing of 
the spheres in a suspension, The viscosity of a suspension is a function 
of the spatial arrangement of the particles in the suspension. Thus any~ 
thing that affects the spatial arrangement would tend to affect the 


viscosity. 


Numerous workers including White and Walton (56), Hudson (25), 
and Furnas (20) have considered the packing of spheres, It is found that 
a bed of m spheres of equal diameter d in the closest (rhombohedral) 
packing contains 74 volume percent solids. In this bed are inter- 
stices of two kinds -''the "square hole" (a reversed spherical cube 
bounded by six spherical surfaces) and the "triangular hole" (a 
reversed spherical tetrahedron bounded by four spherical surfaces). 


Into the square interstices, m spheres of diameter 0.414d can be fitted 
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raising the solid content to about 79 volume percent. Into the triangular 
interstices 2m spheres of diameter 0.225d can be placed at which time 

the solid content is 81 volume percent. This process can be continued 
indefinitely but it can be seen to obtain an arrangement of spheres 
having a high solids content, it is necessary to add spheres of 


different sizes in definite quantities. 


The reduction in the viscosity of suspensions upon the 
addition of fine particles to large particles is then a result of 
the packing of the small particles into the interstices between the 
large particles, which although it increases the volume concentration, 
does not reduce the volume of ''free'' liquid which determines the flow 
behavior of the suspension. Thus a bimodal suspension in which the 
smaller spheres pack into the interstices between the larger spheres 
would be expectéd to have a lower viscosity at the same solid content 
or conversely, would have the same viscosity at a higher total solids 
content. For this reason, it is to be expected that unimodal 
suspensions would tend to have a fairly constant viscosity - 
concentration relationship. Unless it had an extremely wide size dist- 
ribution, a unimodal suspension would have only a few small particles 
to fit into the interstices between the larger particles. There would 
be some effect, most noticeable at high concentrations, but suspensions 
consisting of spheres possessing a bimodal or trimpdal discontinuous 
size-distribution curve would be more likely to SiGe a large reduction 


in viscosity compared to a unimodal suspension. 
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These results may have practical significance in the transport 
of solids in conduits if similar results are found for suspensions of 
materials commonly transported in slurry form. It can be seen from 
Figure 24 that for a bimodal suspension with a diameter ratio of 8.25, 
the viscosity at 40 volume percent is only 65% of the unimodal viscosity 
and at 50 volume percent it is only 36% of the unimodal viscosity for 
the optimum percent of large and small materials. If similar results 
are obtained for practical suspensions, it might be possible to 
prepare the solid phase by segregation of the feed to a slurry 
pipeline to obtain the lowest possible consistency. This will reduce 
the power requirement since the pressure drop is directly proportional 
to the viscosity in the laminar flow region in which a slurry transport 
line might operate. The beneficial effect of reduced viscosity would 
not be as important in the turbulent region because of the reduced 


dependence of the pressure drop upon the viscosity. 


Evidence that this method is being used in industry is given 
in the paper by Dekker and Zimmerman (5), who reveal that in the pro- 
duction of solid rocket fuel, mixtures of various sizes of oxidizers are 
used to lower the consistency of the fuel before casting. It is found 
that the maximum amount of oxidizer is desirable for best performance. 
However, the final mixture must be castable into the varied shapes 
required for rocket motors, It is probable that the work of Sweeny and 
Geckler (45) and Sweeny (46) was done to investigate the possiblility of 


increasing the amount of solid in the slurry while still retaining a low 
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enough consistency so that casting into molds was possible. Dekker and 
Zimmerman show that a mixture of fine and coarse oxidizer in propellant 
in the ratio of 70% coarse and 30% fine particles produces a minimum 
apparent viscosity at the shear rate at which the suspensions were 


tested. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


From this study on the rheology of suspensions of spherical 


particles, it can be stated that: 


i) Over the range of concentration (0 to 50 volume percent) 
and shear rate (0 to 1019 sec” l) Studied the suspensions behaved as 


Newtonian fluids. 


ii) Over the size range studied (47 to 390 microns diameter) 
there is no effect of particle size upon the viscosity of the suspension, 
provided the particles in the suspension are of approximately the same 


size. 


iii) The viscosity of a suspension of spheres is strongly dep- 
endent upon the size distribution of the particles. For bimodal 
mixtures, the minimum viscosity occurs when 70% by volume of the 


total solid phase is composed of large particles. 


iv) Increasing the diameter ratio (the disparity between 
the diameters of the two sizes of spheres in the bimodal mixture) 
decreases the viscosity at concentrations above 15 volume percent. 
The decrease in viscosity is more pronounced the higher the volume 


concentration of solids. 


v) The effect of the size distribution upon the viscosity 


of the suspension can be explained in terms of the packing of the 
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smaller spheres into the interstices between the larger spheres. An 
even greater reduction in viscosity is possible with a selected trimodal 


size distribution of spheres. 


vi) The coefficients in a power series in concentration of 
relative fluidity give a fair correlation with diameter ratio at 


any given percentage of large particles in the suspension. 


vii) The results of this work agree well with those of 
Sweeny (46) who studied the effect of bimodal size distribution 


at volume concentrations of 40 and 55 volume percent solids. 


As an outcome of this work, the following recommendations 
can be made: 

i) A rotational viscometer suitable for the study of the 
rheology of various suspensions should be constructed. A design 
similar to that of the viscometer used by Dudgeon and Michie (9) 
in their work on the rheology of uranium oxide suspensions would 


probably be satisfactory for most suspensions of interest. 


ii) With such a viscometer, a study could be made of the 
effect of size distribution upon the rheology of more practical 
suspensions than spheres. The effect of particle shape and density 
difference between the solid and the suspending solution upon the 


rheological behavior could also be studied. 
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iii) Investigation of the effect of interparticle forces upon 
the rheological properties of suspensions should be made. Since it is 
known that non=Newtonian behavior in suspensions is the result of such 
forces, the measurement of the zeta potential by such a technique as the 
streaming potential and the correlation of this factor as a function of 
the flow behavior of the suspension might yield an insight into the 


mechanism of flow in particle suspensions. 
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NOMENCLATURE 


a Coefficient in power series 

A Designation of 40/60 mesh spheres 

b Coefficient in power series 

B Designation of 100/200 mesh spheres 

c Concentration, volume fraction solids 

¢ Designation of 200/400 mesh spheres 

d Diameter of large spheres in suspension 
du/ dr Shear rate (polar coordinates) 
(du/dr)R, Shear rate at inner cylinder (bob) 

du/ dx Shear rate (cartesian coordinates) 

D Diameter of small spheres in suspension 
Dave Arithmetic average diameter 

Di Diameter of the ith sphere 

1 Mean surface diameter 

D., Mean volume diameter 

DS Mean volume=surface diameter 

e Coefficient in power series 

h Height of cylinder immersed in sample 

k Self-crowding factor in Equation (21) 
ky and ky Instrument constants in Equation (10) for calculation of 


shear rates from speed and dimensions of a viscometer 
K Fluid consistency indéx, defined by Equation (3) 


Ky Shearing stress constant, 4 i= K}® 
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Rate of shear constant for Newtonian fluid, du/dr = KN 


Spring constant of viscometer torsion spring 


Constant 


Torque exerted on the inner cylinder of a concentric cyl- 


inder viscometer 

Flow behavior index, defined by Equation (3) 

Slope of a logarithmic plot of M versus N 

Rotational speed 

Radial distance 

Radius of inner cylinder of a concentric cylinder 
viscometer 

Radius of outer cylinder of a concentric cylinder 
viscometer 

Radius ratio, R,/Ri 

Gap between cylinders of a concentric cylinder 
viscometer 

Linear velocity of outer cylinder of a concentric 
cylinder viscometer 


Mass fraction of material of diameter Ds 


Viscosity of Newtonian fluid 
Viscosity at zero rate of shear 
Relative viscosity = Nees 
Viscosity of suspension 


Viscosity at infinite rate of shear 
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Dial deflection of bob of a: concentric cylinder viscometer 
Crowding factor in Equation (20) 

Coefficient of rigidity of a Bingham plastic fluid 
Kinematic viscosity 

301415 « « « 

Shearing stress 

Shearing stress at bob 

Shearing stress at cup 

Shear rate = du/dr 

Shearing stress at radius r 

Shearing stress at bob 

Yield value of a Bingham plastic fluid 

Relative fluidity = 1WV\ x 

Angular velocity at radius r 

Angular velocity of outer cylinder of concentric 


cylinder viscometer 
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APPENDIX A 


CALIBRATION OF THE VISCOMETER 


Three different cups and three different springs were used 


in this work. 


Bob diameter = aR. = 1,357 inches 
Cup diameter = 2R, 
No. 1 = 1.449 inches 
NoseZz: = 1,507 inches 
No, 3°" = 1.650 inches 
Ss) = Ro /Ry 
Now l= ¥12068 
Nos Zeer) 151105 
Noyes = Sar 216 
Spring Constant = K, 
No. 1 = 387 dyne-cm/degree 
No. 2 = 193.5 dyne-cm/degree 
No. 3 = 77.4 dyne=cm/degree 


For a Newtonian fluid the rate of shear at the bob is: 


(du/dr)p , = 4UN = KoN 
L isis 
The values of Kg are; No. 1 =F O19 
No. 2 = 66.5 
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The shear stress at the bob is given by: 


ae = K,0 = K,@ 
2m Ri oh 


By lumping all the values except 9 into Kj, the end effect is 
absorbed and the value of K, can be determined using standard oils. 


Since for a Newtonian fluid, 


at 


nh (du/er), 


Therefore, 


A 
— 
© 
I 


"KON 
WKoN/ re) 


Ry 


Values of K, determined with four standard oils are given in 
the following table. There appears to be an effect of viscosity and 
gap between cup and bob upon the end effect. The calibration data 


obtained with the standard oils follow the calculated values of Kj. 
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CALIBRATION VALUES OF K) 


he 


shear stress, dynes/cm* 


dial deflection, degrees 


K_9 


il 


OIL VISCOSITY 


1951 cp 45627 


- 522 
2-62 2=62 
2293 Zia DY) 
2.48 2045 


ep 93.6 cp 


Ky 


147.9 cp 


Average Ky 


. 

0.5 

t2.S ca.§ 
: : 


‘ ‘< -- | i : : b a. vl 
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CALIBRATION OF SPRINGS 


SPRING 
NUMBER 


One 


Two 


OIL 


cp. 


56.7 


93.6 


147.9 


TOR 


56.7 


93.6 


147.9 
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Q 
Deflection 


Cup No. 2 


Cup No. 3 


ai: 
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CALIBRATION OF SPRINGS (cont.) 


SPRING OIL N 
NUMBER Cp. Speed 
Three 19,1 10.0 
50) 
BS ey: 
1.67 
5601 10.0 
D320 
3033 
IE(y7/ 
93.46 10.0 
5.0 
Deo 
1.67 
Let 
0.55 


9 
Deflection 
Cup No. 2 Cup No. 3 
128 15 =) 
64 38 
42.5 25 
2055 1245 
- 224 
187 TA, 
2555 7425 
63 37%) 
~ 185 
20245 123 
LOD,.5 6145 
68 41 
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APPENDIX B 


CRITERIA FOR LAMINAR FLOW IN THE ANNULUS OF THE VISCOMETER 


I. Correlation of G.I. Taylor 


CUP Ro C t/Ro logjgt/Ry logio(E)crit (Steric Norit 
NO Vv f AY, 
1 0.06038 0.0038 0.0635 2.8028 37.50 3160 23.4v! 
2 0.06278 0.0063 0.0995 2.9978 3.86 7250 31.7¥! 
3 0.06875 0.0122 0.1776 1.2495 4°25 17806 26.3y! 
Rorrrs radius of outer cylinder, feet 
t = annular gap, feet 
Ve kinematic viscosity, centistokes 
v = kinematic viscosity, eeet ieceend 
U = 2TTNRQ = linear velocity, feet/second 
N = revolutions/second 


II. Gorrelation of H. Schlichting 


ag t/R, WRo)erit NM prantte 
1 0.0635 5 x 10; 23.60! 
2 0.0995 toe 30.5y! 
3 0.1776 se 363i)! 


These correlations are plotted in Figure 30. 
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FIGURE 30 


STABILITY CORRELATIONS 
of 
TAYLOR AND SCHLICHTING 


SCHLICHTING 


Consistent units - ordinates and abscissa dimensionless. 


Oo A OE tS RR RE 


Hert 


hy 


wedi oven 


bine 


O€ SayorT 

CMOTTAINAROD YILITeAT2 
to ; 

MITHOLINOZ CHA POIYAT 


soteelbve - aginn 3oodelaned 


-107- 


APPENDIX C 


CALIBRATION OF THERMOMETER 


Thermometer = Mercury~in-Glass No. 249 

Range: = 1&)%to + 51° € 
Immer sion: 76 mm (Vertical) 

Observed Temperature Standard Resistance 

Thermometer No. 249 Thermometer 

ge ae 

20.00 19.897 
21.00 20.906 
22.00 21.930 
23.03 22.940 
24.00 23.856 
25411 25.005 


These data are plotted on Figure 31. 
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FIGURE 31 


CALIBRATION OF MERCURY=IN=GLASS 
THERMOMETER NO, 249 


RANGE: = 1° to + 51°C 
IMMERSION: 76 mm (VERTICAL) 
NITROGEN FILLED 


21 22 23 


OBSERVED TEMPERATURE 


°o 


C 


24 


25 


26 
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APPENDIX D 


SIZE ANALYSES OF STYRENE-~DIVINYLBENZENE POWDERS 


Material Mesh Size Microns Mass Percent Average Size 
Microns 
A 40/60 500+ 134 545 
pose 420/500 34.3 460 
350/420 20.2 385 
297/350 22.2 324 
250/297 9.5 274 
-250 0.4 230 
100.0 

B 100/200 177+ ney 193 
149/177 bab) 163 

125/149 6.5 137 
105/125 a. 115 

88/105 33.6 96 

74/88 2150 81 
62/74 14.3 68 
53/62 ley 58 

-62 G7 - 

100.0 

Cc 200/400 88+ ey \96 
74/88 lig 'B1 

62/74 153 68 

53/62 18.6 58 

44/53 25.4 48 

37/44 Bed 40.5 
=-37 2 7 ~ 


7110- 


AVERAGE DIAMETERS OF POWDERS 


1. Arithmetic Average Diameter 


Dive = Zs N; D; 


= 2B (x;/D 
& (x; /D;3) 
where Xe = mass fraction of material of diameter Dy 
N; = number of particles of diameter Ds 


2. Mean Surface Diameter 


DNidi2 
| SNi 


D 
Ss 


3. Mean Volume Diameter 


Dd, = 


4. Volume=Surface Diameter 


3 
D = 2,N4Dj = ys (xz) 
ie DNs Wy Gey /D4) 


Using these formulae, the average diameters of the three powders were 


calculated. 


Material 


Mesh Size 


40/60 
100/200 


200/400 


Dave 


356 
1[e) 


38.0 


ue D., 
Microns 
364 BTL 
76.5 7928 
40.1 42.2 


ne 
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APPENDIX E 


EXPERIMENTAL DATA 


The experimental data are presented in the following table. 


The tests are numbered in order, with the letter designating the 


suspension material used. For example, 


A - 40/60 mesh material 
B - 100/200 mesh material 
C - 200/400 mesh material 


BW - 100/200 mesh material after Isopropanol Wash 
Combinations of these designations have been used for the bimodal 
suspensions. As well, the percentages of each material are indicated 
before each series for the bimodal suspensions. In the last two 
series of tests, Numbers 23 and 24, the suspending solution Number 2 
was used. All other tests were conducted with suspending solution 


Number 1. 
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APPENDIX F 


CALCULATED DATA 


The viscosity data of Appendix E are given here as 
relative viscosity and relative fluidity. These were obtained by 
division of the observed viscosity by the viscosity of the suspending 
solution and by division of the viscosity of the suspending solution 


by the observed viscosity respectively. 


Column (4) contains the corrected viscosity. It was 
noted that the observed viscosity of the suspending solution was 
sometimes different from that obtained in the capillary viscometer. 
Since this error was probably due to instrument errors, the observed 
viscosity of all tests in the series was corrected by the addition or 
subtraction of the amount necessary to make the zero percent solids 


value equal to that obtained in the capillary viscometer. 
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(1) (2) (3) (4) (5) (6) 
TEST VOLUME VISCOSITY CORRECTED RELATIVE RELATIVE 
NUMBER PERCENT VISCOSITY VISCOSITY FLUIDITY 
SOLIDS centipoise centipoise Ny 0, 
1=1B 0 1242 = 4 aN - 
1-2B 3438 13.4 = 1.105 0.2906 
1L=3B 508 1449 = lu23 0.815 
1-4B 849 L653 = 1.34 0.745 
1-5B Tlie 3 18.0 - 1.48 0.675 
1-6B 14,2 19<8 = 1.63 0.615 
1-78 1849 23.6 * 1.94 0.4515 
1+8B 2500 3145 - 2.60 0.385 
2=1B 0 12425 7 = ~ 
2=2B 726 15.4 = 127 0.790 
2+ 3B 10.49 17 «2 ” 1.415 0.706 
2-43 14.8 19.8 - 1.63 0.614 
2-5B 18.8 2349 = 1.97 0.509 
2=6B Zouk 274) = 2-26 0.442 
2-78 26.9 3563 ~ 2-90 0.344 
2-8B 3045 45.0 = 3071 0.270 
2+9B 3045 46.4 = 3282 0.262 
2-10B 3461 61.4 _ 5405 0.198 
2=11B 36.4 T1«)d - 6.38 0.157 
2=12B 3845 93.46 - 7e71 0.130 
2~13B 40.5 116 - 9.56 0.1045 
3-1C 0 1343 1242 ae. 7 
3-2C 5x0 15.25 14.15 1. 665 0.858 
3=3C 82/7 171 16.0 1«32 0.759 
3-4C 11.3 18.6 17.5 1.44 0.694 
4-1C 0 1248 Lee ~ ms 
4=2C 4.06 14.4 13.8 1.135 0.880 
4=-3¢ 6.43 154 14.8 1422 0.820 
4-4C 10.2 17«45 16.85 139 0.721 
4=-5C 13.4 1975 19415 LuahS 0.635 
4=6C 16.5 220 Zw 1.76 0.568 
4-7C 19 «7 249 24.3 2200 0.500 
4=8C Z1e/7 2861 2] 5 2426 0.442 
4+9C 25.46 3345 3249 20/1 0.369 
4-10C 21 «2 35 «4 34.8 220/ 0.2349 
4-11¢C 29.45 40.5 39.9 3229 0.304 
4=-12C 3344 5149 5143 4.22 0.237 
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APPENDIX G 


DERIVATION OF EQUATION FOR CONCENTRIC CYLINDER VISCOMETER 


Consider a time~independent fluid of arbitrary behavior 
£(7T) = (du/dr) to be contained in the annulus. Since there is no 
slip at the wall, the liquid in contact with the outer rotating 
cylinder has an angular velocity \%, while that in contact with 
the inner cylinder will be stationary. If the angular velocity is 
uw) at any intermediate radius, the angular velocity gradient is 
dw/dr and the linear velocity gradient across a small element of the 


fluid at the circumference of this circle of radius r is rdw /dr. 


Then rdw/dr = £(T) 
But T(r) = M/(2q r7h) 
rdw/dr =+ (M/ (2th) ) 
At ey Ry; bw = 0 
r= Re wrt 2 


Integrating the above equation between the limits 


a@): = {eosrcon <7h)) dr/x 
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Substituting T(r) = M/ (2? rh) and dr/r = ~dT/(2T) 
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